Non-Abelian Quantum Hall Effect in Topological Flat Bands 
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Inspired by recent theoretical discovery of robust fractional topological phases without a magnetic 
field, we search for the non-Abelian quantum Hall effect (NA-QHE) in lattice models with topological 
flat bands (TFBs). Through extensive numerical studies on the Haldane model with three-body 
hard-core bosons loaded into a TFB, we find convincing numerical evidence of a stable = 1 
bosonic NA-QHE, with the characteristic three-fold quasi-degeneracy of ground states on a torus, 
a quantized Chern number, and a robust spectrum gap. Moreover, the spectrum for two-quasihole 
states also shows a finite energy gap, with the number of states in the lower energy sector satisfying 
the same counting rule as the Moore-Read Pfaffian state. 

PACS numbers; 73.43.Cd, 05.30.Jp, 71.10.Fd, 37.10.Jk 
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Introduction. — The fractional quantum Hall effect 
(FQHE) is one of the most fascinating phenomena in in- 
teracting quantum many-particle systems. While many 
Abelian quantum Hall states have been discovered in a 
two-dimensional electron gas, much effort has been de- 
voted to studies of the non-Abelian quantum Hall effect 
(NA-QHE) since it was proposed two decades ago [ll-[lot. 
One promising experimental candidate for the NA-QHE 
is the 1/ — 5/2 state with electrons occupying the sec- 
ond Landau level (LL) [ll|, while the nature of other 
candidates like the v = 12/5 state remains less settled. 
In addition, the NA-QHE is believed to be possible in 
fast rotating Bose-Einstein condensates (BEG) [l^ ITpl 
and optical lattices with an artificial gauge field [14| . 
The V = 1 bosonic Pfaffian state found in the lowest 
LL [I2I I l3| is favored by a three-body repulsive inter- 
action [3]. In a lattice model, such repulsive interaction 
can be realized by imposing the three-body hard-core bo- 
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son constraint [ij] that no more than two bosons on any 
site are allowed. Interestingly, the three-body hard-core 
bosons can be mapped to spin S = 1 systems, and a 
NA chir al sp in liquid wave function has been proposed 
recently [15|]. The NA-QHE not only is an interesting 
many-body phenomenon, but also provides a promising 
platform for implementing topological quantum compu- 
tation [1^. 

Recently, systematic numerical works demonstrated 
convincing evidence of the Abelian FQHE of interacting 
fermions and hard-core bosons [l7l - [l9j in topological flat 
band (TFB) models Such TFB models belong to the 
topological class of the well-known Haldane model [2lj 
with at least one topologically nontrivial nearly-flat band 
carrying a nonzero Ghern number, which is also sepa- 
rated from the other bands by large gaps 20, 22-24]. 
This intriguing fractionalization effect in TFBs without 
LLs has stimulated a lot of recent research activities 25- 




FIG. 1: (color online). Intensity plots of spectrum gaps in 
U-V phase space for the 20-site lattice at v — 1. E\, E2, 
E3 and Ei denote the energies of the lowest four eigenstates. 
NA, QHE* and SS label the rough phase regions of NA-QHE, 
(non-degenerate) quantum Hall phase, and supersolid inferred 
from all three spectrum-gap plots and other information (see 
the text). 



31|. The V = 1/2 bosonic FQHE found in TFB mod- 
els for hard-core bosons [18] can also be considered as 
one example of the long-sought chiral spin states for spin 
1/2 system [s^- It is now tempting to take the TFB as 
a promising testbed to search for more exotic quantum 
Hall states, possibly with non-Abelian nature. 

In this letter, we search for the non-Abelian quan- 
tum phase for bosons in TFBs without LLs. Through 
extensive exact diagonalization (ED) study on the Hal- 
dane model with three-body hard-core bosons loaded 
into a TFB, we find convincing numerical evidence of 
the 1^ = 1 bosonic NA-QHE, with the characteristic 
three-fold quasi-degeneracy of ground states (GSs) on 
a torus [3, H, [12], an integer quantized Ghern number 
associated with GSs, and a robust spectrum gap in a fi- 
nite region of the parameter space. An energy gap is 
also found to separate the low energy quasihole states 
from the higher energy ones, indicating the existence 
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of the "zero-energy" sector (for the interacting Hamil- 
tonian) [1, [s^, as in the Moore- Read Pfafhan state. The 
number of quasihole states in the lower energy sector 
also satisfy the same counting rule as the Moore-Read 
state 0, H, m, 34 1 . We further obtain the quantum phase 
diagram based on our ED studies and illustrate quantum 
phase transitions of the NA-QHE phase to other compet- 
ing states. 

Formulation. — We study the Haldane model [21] on 
the honeycomb lattice with interacting bosons loaded 
into a TFB E^: 
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where 6j creates a three-body hard-core boson at site r 
satisfying {blf = and {brf = [H; C/ and 1/ are 
the two-body on-site and nearest-neighbor (NN) inter- 
actions. This model can also be considered as a spin-1 
model via the standard mapping from the three-body 
hard-core bosons to the S = 1 spins. It is clear that 
U /t ^ oo corresponds to the limit of the (two-body) 
hard-core bosons. 

The honeycomb lattice has a unit cell of two sites, and 
thus has two single-particle bands. Here, we adopt the 
previous parameters t — 1, t' = 0.60, t" = —0.58 and 
(/)/27r = 0.2, such that a lower TFB is formed with a flat- 
ness ratio of about 50 [3] ■ For our ED study, we consider 
a finite system of Ni x N2 unit cells (total number of sites 
Ns = 2N1N2 and total number of single-particle orbitals 
-^oib = N1N2 in each band) with basis vectors ai and 3.2 
and periodic boundary conditions. We denote the boson 
numbers as N},, and the filling factor of the TFB is thus 

— Nb/Noih- We diagonalize the system Hamiltonian 
in each momentum q = (2Trki / Ni , 2TTk2 / N2) sector, with 
(^1,^2) as integer quantum numbers. 

The phase diagram. — We first look at the spectrum 
gaps for a finite lattice with Ng = 20 sites at filling 1^=1 
as shown by Fig.[Tl where Ei, E2, E3 and E4 denote the 
energies of the lowest four eigenstates, respectively. From 
the three spectrum gaps E4 — E3, E2 — E1 and E3 — E2, we 
can obtain rather rich information of the possible phases 
and related phase diagram. For the 1/ = 1 NA-QHE, 
two necessary conditions are: a ground state manifold 
(GSM) with three quasi-degenerate lowest eigenstates 
(i?3 — El ^ 0), and it is separated from the higher eigen- 
states by a finite spectrum gap E4 — E3 » E3 — E1. From 
the Fig. [ll it can be seen that both conditions are satis- 
fied simultaneously around the left bottom corner in the 
U-V space. The right bottom region is characterized by 
a finite E2 — Ei gap but a very small E3 — E2, which is 
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FIG. 2: (color online), (a) Low energy spectrum E„ — Ei 
versus the momentum feiA^2 -I- ^2 of the NA-QHE phase for 
three lattice sizes Ns = 20, 24 and 28. (b) Scaling plot of the 
spectrum gaps in (a). 



a possible QHE phase (labeled as QHE* with more dis- 
cussion to follow) . While in the upper region with larger 
V, the energy difference E2 — Ei almost vanishes and a 
large E3 — E2 gap appears, indicating the two-fold quasi- 
degenerate states. These are consistent with a sub-lattice 
solid order. Moreover, upon changing boundary phases, 
the two lowest energy states evolve into the higher energy 
spectrum, demonstrating its "metallic" nature in addi- 
tional to its "solid" feature. So we identify this phase 
as a supersolid (SS) phase [35[- We have also obtained 
similar results from a larger lattice with Ns = 24 sites. 

Low energy spectrum. — For the NA-QHE phase, we 
would like to check whether the spectrum gap E4 — E3 
holds for other lattice sizes. A few lowest states in each 
momentum sector of three system sizes with Ng = 20, 24 
and 28 for the case oiU — V — 0.0 are shown in Fig.[5Ja). 
The Hilbert subspaces of the Ng = 28 lattice have the di- 
mensions of about 700 million (which is about the limit of 
the current ED method). We can see that, for each sys- 
tem size, there is an obvious GSM with three-fold quasi- 
degenerate states [two of them in the {ki,k2) = (0,0) 
sector with very close energies]. The GSM is well sep- 
arated from the higher energy spectrum by a large gap 
for all system sizes while the scaling plot of the spectrum 
gap [Fig. [2jb)] suggests that the gap E4 — E3 and the 
three-fold quasi-degenerate GSM of the NA-QHE phase 
should survive in the thermodynamic limit. 

Berry curvature and Chern number. — The Chern num- 
ber [36| (i.e. the Berry phase in units of 27r) of a 
many-body state is an integral invariant in the bound- 
ary phase space [13, HI]: C = ^ J d0id02F{ei,02), 



where two boundary phases 9i and 6*2 are introduced for 
the generalized boundary conditions in ai and a2 direc- 
tions, and the Berry curvature is given by F{6i,92) = 



Im 



\\d92 
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For the GSM of NA-QHE 

phase with Ng — 24, the three GSs maintain their quasi- 
degeneracy and are well separated from the other low- 
energy excitation spectrum when tuning the boundary 
phases, indicating the robustness of the NA-QHE phase 
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(a) U=0.0, V=0.0, NA (b) U=1 .6, V=0.8, QHE* (c) U=1 .6, V=3.0, SS 
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FIG. 3: (color online). Low energy spectra versus 9i at a fixed 
02 = for three phases in the Ns = 24 lattice at v = 1: (a) 
the NA-QHE phase; (b) the QHE* phase; (c) the SS phase. 



(a) U=V=0, E, , Eg (b) U=V=0, (c) U=1 .6, V=0.8, E, 




FIG. 4: (color online). Berry curvatures F(0i, 6*2) A6liA6'2/27r 
at 10 X 10 mesh points for the Ns = 20 cases: (a) the 1st and 
3rd GSs [(0,0) sector] in the GSM of the NA-QHE phase; (b) 
the 2nd GS [(1,0) sector] in the GSM of the NA-QHE phase; 
(c) the single GS of the QHE* phase. 

[Fig. [3Ia)]. Moreover, the GSM in the NA-QHE phase 
shares a total Chern number C = 3: e.g. for the Ns — 20 
cases, we have tvi^o GSs of the GSM in the (0,0) sector 
which contribute the integral Berry phase in [Fig.|3^a)], 
the other GS of the GSM in the (1,0) sector which con- 
tributes the integral Berry phase 27r [Fig.|4l^b)], and thus 
the total Chern number of the GSM is C = 3. In the 
possible QHE* phase, the single GS is well separated 
from other low-energy excitation spectrum when tuning 
the boundary phases [Fig. |31[b)]. The Berry curvature 
of the Ns = 20 case is shown in Fig. HJ^c), which gives 
rise to a quantized Chern number C — 1. On the other 
hand, for the SS phase, the initial two-fold GS quasi- 
degeneracy is immediately destroyed when tuning the 
boundary phases, the two GSs evolve into the higher ex- 
citation spectrum. They do not have well-defined Chern 
numbers since there is no well defined spectrum gap, 
which indicates a "metallic" feature of this SS phase 
[Fig. He)]. 

Quasihole excitation spectrum. — In order to investigate 
the possible fractional statistics of the NA-QHE state, we 
study the quasihole excitation spectrum by removing one 
boson from the v = 1, and expect two quasiholes of frac- 
tional bosonic charge 1/2 [l|,[l,|3]- As shown in Fig.[5l[a), 
for a typical NA-QHE state on the Ns — 24 lattice, the 
quasihole spectrum exhibits a distinguishable gap which 
separates a few lowest states in each momentum sector 



(a) U=0.0, V=0.2, NA (b) (ki,k2)=(1,1) (c) (ki,k2H2,0) 
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FIG. 5: (color online), (a) Quasihole excitation spectrum in 
the NA-QHE phase for Ns = 24 and iVb = 11. (b)-(c) Low 
energy spectra versus 9i at a fixed 02 = in two momentum 
sectors of (a). 
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FIG. 6: (color online). Root configurations in the A'orb = 12 
single-particle orbitals: (a)-(c) three GSs (02), (20) and (11); 
(d)-(f) two-quasihole states with two domain walls (repre- 
sented by two vertical lines). 



from the other higher-energy states. For each momen- 
tum sector, we check the spectrum gap upon changing 
boundary phases. For the sectors already with a large 
gap [e.g. the (1,1) sector], the spectrum gap is main- 
tained well for all boundary phases with 6 lowest states 
below the gap as shown in Fig. [SJb) . For other sectors 
[e.g. the (2, 0) sector] where the quasihole gap is less ob- 
vious, upon changing boundary phases the spectrum gap 
becomes clearer and we also find 6 lowest states below 
the gap as shown in Fig. \^c). By summing up all 12 
sectors together, we have 72 low energy quasihole states 
in total. We also find similar features for the Ns = 20 
(and Nb = 9) case: there are 5 quasihole states in each 
momentum sector, and all 10 sectors give 50 low energy 
quasihole states in total. 

The number of low energy two-quasihole states in the 
NA-QHE phase described above can be heuristically un- 
derstood from the counting rule based on the general- 
ized Pauli principle M. |l9l]. Using the Wannier repre- 
sentation for a TFB '25|, a set of iVorb = A^s/2 peri- 
odic single-particle orbitals are formed. Now, let us con- 
sider a system with iVorb = 12 as an example. The gen- 
eralized Pauli principle that no more than two bosons 
occupying two consecutive orbitals 0, @, [l3| results in 
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(a) E vs. U, at V=0.0 (b) E vs. V, at U=0.0 (c) E vs. V, at U=1 .6 




U V V 



FIG. 7: (color online). Quantum phase transitions when tun- 
ing U (or V) with V (or U) fixed for N, = 20: (a) NA-QHE 
to QHE*; (b) NA-QHE to SS; (c) QHE* to SS. 



only three GS root configurations in the above orbitals 
\nx„nx,,...,nx^ J: (02) = |020202020202), (20) = 
1202020202020), and (11) = |111111111111) [Figs. IJa)- 
(c)]. We now count how many ways we can re- 
move one boson from the above three GS configura- 
tions (02), (20) and (11). The boson occupancy of 
two-quasihole states should be a mixture of two seg- 
ments of the three GS configurations, with two do- 
main walls each representing one quasihole with 1/2 
charge Q- A simple analysis gives 6 types of con- 
figurations with odd number of I's: |. . . 20|1|020 . . . ), 

|...20|111|020...), I ...20|11111|020...), , and 

|0|11111111111|) [Figs.|6i;d)-(f)], where two domain walls 
(quasiholes) are displayed by two vertical lines |'s. Con- 
sidering 12 translations of the above 6 states, we finally 
get the total 72 in general) two-quasihole states 

in exact accordance with our numerical results. 

We would like to emphasize that, for the GS obtained 
numerically, the above mentioned root configurations are 
not the dominant configurations due to quantum fluctu- 
ations just like the other FQHE states in LLs on a torus. 
But the number of the low energy sector quasihole states 
in QHE systems is a consequence of the Pauli princi- 
ple [tI IsI, u3\ resulting from strong short-range pseudopo- 
tentials j9„ ^3]. Indeed, based on the pattern of zeros 
classification for FQHE states, the Moore- Read state has 
a lower energy sector with the same number of quasiholes 
as obtained in our numerical results {sj]. 

Quantum phase transitions. — When tuning U (or V) 
with fixed V (or U) away from the NA-QHE region, we 
observe quantum phase transitions from the NA-QHE to 
other quantum phases including a possible QHE* and the 
SS as shown by Fig. [71 

Even though the NA-QHE phase is shown to be re- 
markably robust in our ED study, we are less certain 
about the QHE* phase. As shown in Fig. [7l[a), if we in- 
crease [/, two higher energy states from the GSM of the 
NA-QHE phase will emerge into the excited spectrum, 
while the evolution of the wave function of the lowest 
energy state is very smooth. Thus we suspect that the 



system may have a very long correlation length near the 
possible QHE* phase, and we conjecture that the spec- 
trum gap above the GS may collapse when the system 
size becomes very large. We leave this issue to be ad- 
dressed in future studies. 

The SS phase has a two-fold GS quasi-degeneracy, 
strong intra-sublattice density-density correlations and 
vanishing inter-sublattice density-density correlations. 
These observations indicate that the bosons prefer occu- 
pying one of the two sublattices. By changing boundary 
phases, the two GSs of the SS state evolve into the higher 
energy spectrum, indicating its "metallic" feature besides 
its solid feature. However if we fix the V as 3.0 and go to 
a larger U > 8.0 (not shown in the phase diagram), the 
"metallic" character of the SS phase disappears while the 
solid feature remains indicating the system enters a pure 
solid phase. 

Summary and discussion. — The v — 1 NA-QHE is 
characterized by a three-fold quasi-degenerate GSM and 
a total Chern number C = 3 carried by these states. 
The three-fold quasi-degenerate GSM is only observed 
for system with even number of bosons {Ni, — 10, 12, 14) 
a,t u ~ 1, while it is absent for odd numbers of bosons 
(e.g. Nf, = 9 and Ng = 18) consistent with the pairing 
nature for Pfaffian-like states. Interestingly, the ly — 1 
NA-QHE is already quite stable with a clear GSM and a 
large spectrum gap for three-body hard-core bosons with- 
out additional interactions {U = = 0). The spectrum 
gap can even be significantly enhanced with the presence 
of a small V, which makes it possible for such a state to 
be realized using optical lattices. 
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Note added. — After the completion of this work, we 
became aware of a related work by Bernevig and Reg- 
nault [39] addressing the NA-QHE in fermionic TFB sys- 
tems. 
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